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We prove the existence of a generalization of Kelvin's circulation theorem 
in general relativity which is applicable to perfect isentropic magnetohydro- 
dynamic flow. The argument is based on a new version of the Lagrangian for 



, perfect magnetohydrodynamics. We illustrate the new conserved circulation 



with the example of a relativistic magnetohydrodynamic flow possessing three 
symmetries. 



I. INTRODUCTION 



> 

It is hard to overstate the important role that Kelvin's theorem on the conservation 
O ! of circulation of a simple perfect fluid has played in the development of hydrodynamics. 

Among other things it provided the basis for the discussion of potential flows, and showed 
that isolated vortices should exist, and that they obey the Helmholtz laws, etc. On the 
other hand, Kelvin's theorem is fragile: as soon as dissipation comes in, it breaks down. 
And when the body force per unit mass of fluid is not a gradient, as happens for the Lorentz 
force in magnetohydrodynamics (MHD), Kelvin's theorem ceases to apply. 



£jrj! Often fluids in the real world carry magnetic fields. For example, the fluid at the Earth's 

core, the plasma in the sun and pulsars, the ionized gas in interstellar space and in su- 
pernova remnants, the plasma in intergalactic space, and many others, all carry magnetic 
fields. These are important systems for which we need insights of the sort Kelvin's theorem 
bestowed on ordinary fluid dynamics. Can we extend Kelvin's theorem to MHD ? 

The first such extension was found two decades ago by E. Oron fl working with the formal- 
ism of relativistic perfect MHD. This is a circulation theorem for a hybrid velocity-magnetic 
field combination. Oron's derivation assumes both stationary symmetry and axisymmetry, 
while it is well known that Kelvin's theorem requires neither of these. Yet it has not proved 
possible to rid Oron's result of the symmetry assumptions. 

In the present paper we follow, on the wake of our earlier paper, @ a different route. We 
use the least action principle to give a rather straightforward existence proof for a generically 
conserved hybrid velocity-magnetic field circulation in general relativistic MHD which does 
not depend on the presence of spacetime symmetries. Recently Elsasser B has given a related 
result which he obtains by direct manipulation of the relativistic MHD equations. 

A formal introduction to relativistic MHD is given by Lichnerowicz. As mentioned, 
we approach the whole problem not from equations of motion, but from the least action 
principle. In special relativity Penfieldil proposed a perfect fluid Lagrangian which admits 
vortical isentropic unmagnetized flow. The early general relativistic Lagrangian of Taub @ as 
well as the more recent one by Kodama et. al'i describe only irrotational perfect fluid flows. 
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The Lin device B to include vortical flows is incorporated by Schutzo in his Lagrangian. 
Carter has introduced a relativistic Lagrangian for particle-like motions from which the 
properties of fluid flows, including vortical ones, can be inferred. However, it does not 
correspond in detail to the MHD paradigm. AchterbergQ proposed a general relativistic 
MHD action, which, however, describes only "irrotational" flows. Thompson used this 
Lagrangian in the extreme relativistic limit. Heyl and Hernquist © modified it to include 
QED effects. In this paper we follow Schutz'so approach while supplementing it by the 
introduction of magnetic fields 

In Sec. II. A we review the relativistic MHD equations. In Sec. II. B we describe our 
Lagrangian, @ and show that it gives the correct Maxwell and fluid equations, while in 
Sec. II. C we recover the relativistic Euler equation from it. In Sec. III. A we derive the 
general form of the conserved circulation, while in Sec. III. B we illustrate it with the special 
case of a MHD flow endowed with three symmetries. 



II. RELATIVISTIC ACTION PRINCIPLE 

In this section we construct a Lagrangian density for MHD flow in general relativity 
(GR). Greek indices run from to 3. The coordinates are denoted 

x° stands for time. A comma denotes the usual partial derivative; a semicolon covariant 
differentiation. Our signature is {—,+,+, +}. We take c = 1. 



A. Relativistic MHD Equations 

Our first step is enumerating all the correct general relativistic (GR) equations for MHD. 
These were developed by Lichnerowicz, @ Novikov and Thorne, @ Carter, @ Bekenstein and 
Oronfl and others. 

The first equation states the conservation of the number of particles (we do not consider 
particle annihilation or creation processes), 

N a , a = (nu a ) ]a = 0, (2.1) 

where N a is the particle number 4-current density, n the particle proper number density 
and u a the fluid 4-velocity field normalized by u a u a = —1. We consider flows which are 
inviscid and adiabatic, and therefore s, the entropy per particle, is conserved along flow 
lines: 

(sN a ). a = or u a s, a = 0. (2.2) 

The energy momentum tensor for the magnetized fluid is obtained by adding the elec- 
tromagnetic energy-momentum tensor to that of an ideal fluid: 

T ap = pg *p + { p + p ) U <* U P + (4 7 r)- 1 (F a7 F /3 7 - -F^F lS g aP ). (2.3) 

4 

Here p denotes the fluid's proper energy density (including rest mass) and p the scalar pres- 
sure (assumed isotropic), while F al3 denotes the electromagnetic field tensor which satisfies 
the GR Maxwell's equations 
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F a P.p = 4vrj Q (2.4) 

FaPtf + F/3j,a + F-ya,/3 = 0, (2.5) 

where j a denotes the electric 4-current density. 

The magnetic Euler equation for the fluid is derived from the vanishing covariant diver- 
gence law T a/3 .^ = 0: 

(p + p)u p u a , p = -(g a ? + u a u?)p^ + (4vr)- 1 F a/3 F/ ;7 . (2.6) 

The term a a = u^u 01 ^ stands for the fluid's acceleration 4-vector. The use of covariant 
derivatives and curved metric ensures that effects of gravitation are automatically included. 
In view of Eq. ( |2.4|) the above yield the MHD Euler equation 



( p + p)a Q = -^% + F^, (2.7) 
where h af3 is the projection tensor 

The Euler equation describes a general electromagnetic field carrying flow. One needs 
an additional condition to distinguish MHD flow from all others. For any flow carrying an 
electromagnetic field, the (antisymmetric) Faraday tensor F a ^ may be split into electric and 
magnetic vectors with respect to the flow: 

E a = F af} v? (2.9) 
B a = *F l3a u>' = l -ep ai 8 F 7 V. (2.10) 

Here e a ^s is the Levi-Civita totally antisymmetric tensor (eoi23 = (—g) 1 ^ 2 with g denoting 
the determinant of the metric gap) and *F a/ 3 is the dual of F a p. In a frame comoving with 
the fluid, these 4-vectors have only spatial parts which correspond to the usual E and B, 
respectively. One can use Eqs. ( |2.9| - |2.1CI| ) to express F a p using those 4-vectors 

F af3 = u a E p - u p E a + e^su'B 5 (2.11) 

For an infinitely conducting (perfect MHD) fluid, the electric field in the fluid's frame must 
vanish, i.e., 

E a = F a p U P = 0. (2.12) 
This corresponds to the usual MHD condition E + v x B = 0. 

B. Relativistic Lagrangian density and equations of motion 

We now propose a Lagrangian density for GR MHD flow based on Schutz'sil Lagrangian 
density for pure fluids in GR: 

C = -p(n, s) - (IGtt)- 1 ^^ + cf>N a , a + V (sN a ). a + A (jN a ). a + r a F a ^. (2.13) 
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As shown below, our Lagrangian reproduces Eqs. ( |2.1| - |2.2| ), (|2.4[ - |2.7|) and ( |2.12|) . Here <fi is 
the Lagrange multiplier associated with the conservation of particle number Eq. (|2.1|) viewed 
as a constraint, rj is that multiplier associated with the adiabatic flow constraint, Eq. ( j2.2| ), 
and A is that associated with the conservation along the flow of Lin's I qnantity 
a quartet of Lagrange multipliers which enforce the field freezing condition Eq. ( [2.12j ). We 
should interject that r a is not determined uniquely, as we discuss in Sec. III. 

We view 7, N a and s as the independent fluid variables, while n and u a are determined 
by the obvious relations 

- N a N a = n 2 ; u a = n" 1 N a . (2.14) 

Some authors prefer to include in the Lagrangian the constraint N a N a + n 2 = 0, which 
stands for the normalization of the fluid's 4-velocity. However we choose to impose this 
constraint later and thus remain with a simpler Lagrangian. 

We can now vary the Lagrangian with respect to the independent variables. Variation 
of 4> recovers the conservation of particles N a - a = 0. Variation of A with subsequent use of 
the previous result yields 

7 iQ m° = 0. (2.15) 

If we vary Lin's 7 we get 

A, a w Q = 0. (2.16) 

These results inform us that 7 and A are both locally conserved along the flow. In view of the 
thermodynamic relation n' 1 (dp/ds) n = T, with T the locally measured fluid temperature, 
variation of s gives 

u a r], a = -T. (2.17) 

We now vary N a using the obvious consequence of Eq. (|2.14|) , 

5n=-u Q 5N a } (2.18) 

together with the thermodynamic relation @ involving the specific enthalpy /i, 

fi = (dp/dn) s = n~ l (p + p). (2.19) 

We thus get the most important equation herein: 

fJ>u a = <p, a + srj >a + 7A, Q + T p F a p. (2.20) 

By contracting Eq. ( |2.20| ) with u a and using u a u a = — 1 as well as Eqs. ( |2.12| ) and ( p. 16 - 
ET7D , we get 

cj) a u a = -n + Ts. (2.21) 

Thus the proper time rate of change of <fr along the flow is just minus the specific Gibbs 
energy or chemical potential. 
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The importance of using Lin's 7 is clear from Eq. (|2.20|) . In the pure isentropic fluid 
case [F a/3 = and s = const.), the Khalatnikov vorticity tensor given by 

UJ a( 3 = (/^) )Q ~ 0«a) j/3 = {l\p) 7a ~ (l\a) jP (2.22) 

would vanish in the absence of 7, thus constraining us to discuss only irrotational flow. This 
problem is well known from non relativistic pure fluid Lagrangian theory. Lin 11 remarked 
that one can label each fluid element by its original Lagrangian coordinate. The requirement 
that this stay fixed adds an additional constraint ("label conservation") to the Lagrangian 
function, and makes possible the description of vortical flow. While for isentropic flow 
Kelvin's theorem forbids the creation of vorticity, the flow in any given region can be vortical 
due to conditions upstream. 

As customary, we write F a/3 = Ap. a —A a ./3 = A/3 j0l —A a> p, which ensures that the Maxwell 



Eqs. (2J3) are automatically satisfied. The other half, Eqs. (|2.4j ), are obtained by varying 
with respect to the components of the vector potential A a . Because of the antisymmetry 
of F a p, the last term of the Lagrangian, Eq. (|2.13|) , can be written as {r f5 N a — r a N l3 )A a ^. 
The variation of A a in the corresponding term in the action produces, after integration by 
parts, the term [{-g) x / 2 {T a N p - T P N a )}^5A a . Because for any antisymmetric tensor t ap , 
{-g) l/2 t al3 .p = [(-g) 1/2 t Q/3 ] i/3 , variation of A a leads to the equation 

F^.p = 4vr (r a N p - T?N a ) . (2.23) 



We see that this is just Eq. ( |2^) provided we identify the electric current density j a as 

f =(r<*N 13 -r p N a ) (2.24) 

Since the divergence of the divergence of any antisymmetric tensor vanishes, the charge 
conservation equation (j a - a = 0) is satisfied automatically. Formally Eq. (|2.23j ) determines 
the Lagrange multiplier 4- vector r a , modulo the freedom inherent in it, as we discuss in 
Sec. III. 



C. MHD Euler equation in General Relativity 

We now go on to tie the equations of motion together to yield the MHD Euler equation 
(2.7). We begin by writing the Khalatnikov vorticity uop a in two forms, 

Ufia = fi,pU a - fi^Ufs + flU a;l3 - /lUfca, (2.25) 

as well as by means of Eq. ( |2.20| ) 

+ T 5 ;(3 F a5 - T & ]a Fp 5 + T S F a5 .(3 - T S FpS;a- (2.26) 

Contracting the left hand side of the first with N a , recalling Eq. ( |2.14 ) and that by nor- 
malization u a u a]l 3 = whereas vPu a .fi = a a (recall that a a is the fluid's 4-acceleration), we 
get 
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-n\ifi — n/j, ja u a Ufs — niiap = —nhp a p, )a — npiap. 



[2.27) 



On the other hand, contracting Eq. ( 2.26|) with N a and using Eqs. ( p. 15 - 2. 17 ) and (p. 12 ) 
to drop a number of terms, we get 

u Pa N a = -nTs j/3 - T S . a Fp$N a + r 5 F a5 , p N a - r 5 F p5 , a N a . (2.28) 

—nTsfi is the same as —nThn a s n . It is convenient to use the 
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By virtue of Eq. 

thermodynamic identity dp = n _1 dp + Tds, which follows from Eq. ( |2.19| ) and the first law 
d(p/n) = Tds — pd{\/n), to replace —nTs^ in Eq. fl2.28|) by hp a (—np ja +p jQ! ). Equating 
our two expressions for u)p a N a gives, after a cancellation, 



- {rijiap + h/3 a p,, 



-T S , a Fp 5 N a + r 5 F aS .pN a - T S F ps , a N a . 



(2.29) 



The last two terms in this equation can be combined into a single one by virtue of 
Eq. (|2.5|) , which can be written with covariant as well as ordinary derivatives. Further, by 
Eq. ( |2.19D we may replace n/i by p + p. In this manner we get 



(p + p)a p = -h p a p, a + F Pa , 5 T & N a + F p& r d , a N' 



(2.30) 



The term r s . a N a here can be replaced by two others with help of Eq. (|2.23|) if we take into 
account that N 13 .^ = 0: 

(P + P) = -h p a p, a + {A7t)- l Fp S F 5a , a + F Pa , 5 T 5 N a + F ps (r a N s ) ^ . (2.31) 

We note that the last two terms on the right hand side combine into (^Fp a N a T 5 ^j which 

vanishes by Eq. (|2.12 ). Now substituting from the Maxwell equations (|2.4|) we arrive at the 
final equation 



(p + p)a p = -hp a p ja + Fpsj* 



(2.32) 



which is the correct GR MHD Euler equation (|2.7| ). 

Note that we have not used any information about r a beyond Eq. ( ^.23[ ); hence Euler's 
equation is valid for all choices of r a , of which there are many as we shall explain. Since 
we are able to obtain all equations of motion for GR MHD from our Lagrangian density, we 
may regard it as correct, and go on to look at some consequences. 



III. NEW CIRCULATION CONSERVATION LAW 



A. General Remarks 



Eqs. ( j2.20| ) and ( 2.15 - 2.16 ) lead immediately to a law of circulation conservation for 
relativistic perfect isentropic MHD flow. Define the vector field 



Z a — pu a F a £ 



(3.1) 
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and its associated circulation T 



z a dx a , (3.2) 



where C is a closed simply connected curve drifting with the fluid. According to Eq. fl2.20| ) 



z a = 4>,a + sr},a + 7-\ Q . Since (f) :a is a gradient, its contribution to T vanishes. Likewise, for 
isentropic flow (s = const.) the term involving sr], a makes no contribution to V. Thus 

7A a dx a = f 7 dX. (3.3) 
c ' Jc 

By Eqs. ( |2 . 1 5| - pTT6|) both 7 and A are conserved with the flow. Thus T is a circulation which 
is conserved along the flow. 

In the absence of electromagnetic fields and in the nonrelativistic limit (/^ — > m where 
m is a fluid particle's rest mass), T for a curve C taken at constant time reduces to Kelvin's 
circulation. On this ground our result can be considered a proof that a generalization of 
Kelvin's circulation theorem to GR MHD exists. This conclusion goes beyond Oron's original 
conserved circulation in MHD0 in that no symmetry is necessary here for the circulation 
to be conserved. 

To make full use of the new conservation law to solve or simplify problems in MHD, one 
must evidently know r a explicitly. There is a certain amount of freedom in r a which we 
have already discussed, il Here it is important that the Lagrangian density fl2.13|) is invariant 



under the addition of fu a to r a , where / is an arbitrary scalar, because fF a pu a viP vanishes 
identically. We use this freedom to demand that r a u a = 0. We now recast the circulation 
law in a form eminently suitable for use in problems with symmetry. 
First from Eq. ( |2.23| ) we infer that 

T a N P _ T P N a = (4^-1(^/3 _ W*f)} . W af3^ = q (3 ^ 

where W a/3 is an antisymmetric tensor. Since W a ^ must be divergenceless, it can generically 
be written as the dual of a curl, 

where = A$ r/ — Ay t s. A a here is to be distinguished from the ordinary vector potential A a 
of F a p. Parenthetically we mention that one can independently make gauge transformations 
of A a and of A a . This underscores the little known fact that MHD is a theory with U(l) x 
£7(1) gauge symmetry. 

By taking the dual of Eq. ( |3.4|) and then contracting it with u 1 (remembering that twice 
dual is equivalent to changing the sign) we get, with the help of Eq. ( |2.10| ), 



Bs = T^v? . (3.6) 
Contracting Eq. ( |3.4|) with up and recalling that F a ^up = and r a u a = gives 

T a = (^nyh^T^up. (3.7) 
Now using Eq. ( [2.1 1|) with E a = and Eq. (|3.7j ) we have 
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FcprP = (STrnyh^se^u^B 5 ^^ (3.8) 
With the easily checked identity 

-5^ a 5^5 u s + S" a 5\5^ - S u J%^s (3.9) 

Eq. ( |3.8|) reduces to 

F aP r p = {^n)-\T av B v - B v B v u a ) (3.10) 

where use has been made of Eq. (|3.6| ), B^u u = and u^u^ = — 1. The conservation law 
( p.2| ) now takes the form 

T = £ [ X u a - {^n)- 1 T aP B l3 ]dx a 1 (3.11) 

with x = A*+ {4:7m)~~ 1 Bf3B t3 . The quantity x plays an important role in Oron's generalization 
of Kelvin's theorem to MHD. fl 



B. Example: MHD Flow with Three Symmetries 

Since T a p is not known explicitly, we cannot work out the conserved circulation without 
further work. Here we shall make some progress in this direction in the case of flow which 
is both stationary, and possesses two additional spatial symmetries (these last we assume 
not to be of the angular type). This means the physical quantities such as u a , N a , j a 
or B a are unchanged upon being Lie dragged along either of the three Killing vectors. 
We cannot automatically require the same of T a p because it is not a directly measurable 
quantity. However, as mentioned in our earlier work, o we can make a transformation 
J~ap ~^ J~ a j3 + f a /3, where f a p is a curl and orthogonal to u 13 , without changing the values of 
j a or B a ; this transformation at most adds to T a conserved quantity leaving it conserved. 
We shall assume here that by means of such a transformation we can make JF Q(3 share the 
symmetries of j a or B a . 

The following remarks apply when we choose the coordinates x 2 and x 3 to extend along 
the integral curves of the two spatial Killing vectors; by assumption these curves are non- 
compact, and so are x 2 and x 3 . We also assume x 1 is noncompact. The most general form 
for the "vector potential" A a for which T a p is independent of x°, x 2 and x 3 , is (here we 
sacrifice manifest covariance in order to make the symmetries manifest) 

A a = x°<l a + x 2 ^ iQ + x 3 H Q + V a (3. 12) 

where ^ and H are each a linear combination of x°,x 2 and x 3 with constant coefficients 
plus a function of the nontrivial coordinate x 1 only, while the components of the "vector" 
V a also depend only on x 1 . Thus 

Foi = ($ - V ),i = $ i (3.13) 
^21 = $-V 2 ),i = V,i (3.14) 
J r 3 i={E-V 3 ), 1 = E tl (3.15) 
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while ^02, ^03 an d T 23 are all strictly constant, and thus stand for global parameters of the 
flow. 

By means of Eq. ( |3.6| ) we may now compute the components of B a : 

B = $ V + F 02 u 2 + F 03 u 3 (3.16) 

B 1 = - ^ , 2 m 2 - E tl u 3 (3.17) 

B 2 = # X + F 20 u° + F 23 u 3 (3.18) 

B 3 =S/+J 3 o« + ^ 2 (3.19) 



Solving these last for the derivatives of $, \1/ and S we get from ( |3.13| )-( ]3.15| ) 



F m = (B - F Q2 u 2 - F G3 u 3 )/u l (3.20) 
F21 = (B 2 + F 02 u° - F 23 u 3 )/u l (3.21) 
^•31 = (B 3 + F 03 u° + T 23 u 2 )ju x (3.22) 

We now calculate the quantity appearing last in Eq. ( |3.11| ) with the help of Eqs. ( |3.20| )- 
Q3.22j ). If temporarily we take Tq 2 = JF 03 = 0, we get 

F a pB p dx a = (B^u^iBodx + B 2 dx 2 + B 3 dx 3 ) - (B B° + B 2 B 2 + B^B^dx 1 /u 1 ) 

+ (F 23 /u l )[(u :i B 2 - ^B^dx 1 + (^B 3 - u 3 B l )dx 2 + (^B 1 - u l B 2 )dx 3 } (3.23) 

This is much simplified by adding BiB l (dx l /u l ) to the first term and subtracting it from 
the second term. In addition, one can unify the terms in square brackets by employing the 
Levi-Civita tensor. Putting all this together and using the expression ( [2.11| ) we have 

T = £ [ X u a + (Airnu 1 )- 1 (B^5\ - B x B a + (F 23 /^)F 0a )} dx a (3.24) 

We now show that the last term in the integrand does not contribute. First of all we can 
take the constant T 23 out of the integral. Next we realize that the law of particle number 
conservation (|2.1| ) and the assumed symmetries tell us that — C^ 1 = ^/—gnu 1 does not 
depend on x l either, and so can also be taken out of the integral. We are left with a term 
proportional to C §F 0a dx a . Now, just like jF a/3 , F a p derives from a vector potential, and 
must have all the symmetries we have assumed. We can thus take its vector potential A a 



to have the form ( |3.12| ), but with new functions (sans circumflex), each of which 

is, again, a linear combination of x°,x 2 and x 3 with constant coefficients plus a function of 
the nontrivial coordinate x l only. It follows that F 0a dx a = Q'^x 1 + const. dx 2 + const. dx 3 . 
But this is a perfect differential; hence the term proportional to T 23 in Eq. ( [3.24D vanishes. 
If we now reinstate JF 02 and JF 03 we find that they contribute to T terms proportional to 
CFq 2 § F 3a dx a and CJF 03 § F 2a dx a , both of which are found to vanish by reasoning analogous 
to the above. 

By the symmetries the term B pB 13 [nu l )~ l in the integrand of Eq. (|3.24j) can only depend 
on x 1 . It is integrated over x l only, once forward and once backward because C is a closed 
curve. Hence this term makes no contribution to T. Further, it is a consequence of Euler's 
equation ( |2.7| ) that (fxB a ). a = oM Only the radial derivative survives here, so that we have 
^/—gB 1 = D//J, with D a constant, a further global parameter of the flow. We thus obtain 
the final form of the conserved circulation: 
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Xu a + DC (Ait dx a (3.25) 

This is exactly Oron's original conserved circulation @ for two symmetries (we have here 
defined C and D to correspond to the two quantities with the same names which are con- 
served along streamlines for the case of two symmetries ) . Thus the imposition of a third 
symmetry does not cause any changes in V. In future work we shall endeavor to recover 
Oron's conserved quantity in the presence of two symmetries from the present approach, 
as well as to explore the conserved circulation in stationary flow with no spatial symmetry. 
The general case of no symmetries is a more distant goal. 
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